Proposed method for direct measurement of non-Markovian character of the qubits 

coupled to bosonic reservoirs 
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The non-Markovianity is a recently proposed characterization of the non-Markovian behavior in 
an open quantum system, based on which we first present a practical idea for directly measuring 
the non-Markovian character of a single qubit coupled to a zero-temperature bosonic reservoir, 
and then extend to investigate the dynamics of two noninteracting qubits subject to two reservoirs 
respectively with a lower bound of non-Markovianity. Our scheme, with no need of optimization 
procedures and quantum state tomography, is helpful for experimental implementation. 
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Realistic quantum systems are fragile due to unavoid- 
able interaction with the environment [l|, For this 
reason, the dynamics of open quantum systems have 
attracted much attention in the investigations of mod- 
ern quantum theory, particularly of quantum informa- 
tion processing (QIP) Over past decades, the con- 
ventionally employed Markovian approximation with the 
assumption of infinitely short correlation time of the en- 
vironment has experienced more and more challenges 
due to advance of experimental techniques 1] , and non- 
Markovian features have been observed in some physical 
systems, e.g., atomic and molecular systems high-Q 
cavity systems @ and solid state systems @. To under- 
stand these non-Markovian effects, there have been vari- 
ous kinds of analytical and numerical methods developed 
so far [Ul7l-[Toj. such as non-Markovian quantum trajecto- 
ries Q, pseudomodes Q, non-Markovian quantum jumps 
@, and quantum semi-Markov processes [Toj . 

Among the recent investigations for non-Markovian 
behavior [lTl - flij ]. a particular approach using non- 
Markovianity contains the exact characterization of the 
non-Markovian behavior of a quantum process without 
making any approximation for the dynamics [l2j . The 
non-Markovianity is defined to quantify the total amount 
of information flowing from the environment back to the 
system through a quantum process <&(t) with p(t) = 
$>(t)p(0), where p(0) and p(t) are the density operators 
of the system at initial time and at arbitrary time. In 
Ref. [12j, the information is characterized by the trace 
distance V (pi,p 2 ) = \tf \px — P2 1 of two quantum states 
Pi and p 2 , describing the distinguishability between the 
two states, and satisfying < V < 1 [3]. T> maximally 
reaches 1 when the two states are totally distinguishable 
and approaches for two identical states 0]. The di- 
rection of information flow is dependent on the slope of 
V( Pl (t),p 2 (t)) , i.e., when d t V { Pl {t) , p 2 {t)) < 0, the in- 
formation dissipates to the environment and vice versa. 



Therefore, the non-Markovianity could be calculated by 



Pl,2(0) 



E 



ts( /_nuil\ / r 



: ),P2 (rr x )) 

l ))L 



(i) 



with r™ in (r™ ax ) the time point when V ( Pl (t), p 2 (t)) 
reaches the nth local minimum (maximum). Eq. (1) can 
be carried out by summing up the amount of the increase 
of the trace distance over each time interval [r™ ln , r™ ax ] 
for any pair of initial states pi(0) and p 2 (0), where the 
maximum is considered as the non-Markovianity jV($). 
Since this is a problem of optimization, however, we have 
to consider all pairs of initial states in our calculation, 
which is inconvenient and impractical, especially from 
the viewpoint of experimental exploration. 

In this work, we show that the optimization prob- 
lem of Af(&) could be simplified to an effectively com- 
putable expression in the case of a single qubit coupled 
to a zero-temperature bosonic reservoir. Moreover, for 
two independent qubits coupled to two bosonic reser- 
voirs respectively, we investigate a lower bound of the 
non-Markovianity. The favorable feature of our method 
is the possibility to connect the non-Markovianity to the 
population of the qubit in the excited state, which could 
be detected directly in experiments without the require- 
ment of tomographic reconstruction of the density ma- 
trix. In addition, our result also makes it possible to 
have an easy evaluation of non-Markovianity even with- 
out much information about the interaction between the 
qubit and the reservoir. 

We first consider a single qubit coupled to a zero- 
temperature bosonic reservoir. The Hamiltonian in units 
of h = 1 is given by 



H =u Q \e) 
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where ujq is the resonant transition frequency of the qubit 
between the excited state |e) and the ground state \g). 
wi and ai (aj) are, respectively, the frequency and the 
annihilation (creation) operator of the Ith mode of the 
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reservoir with the coupling constant gi to the qubit. The 
dynamics of the single qubit can be represented by the 
reduced density matrix [l| 



A*) 



P?e(0)\Ht)\ 2 P ? g (0)b(t) 

pl(0)b*(t) l-p£(0)|6(t)| s 



(3) 



in the qubit basis {|e) , \g)} , where the superscript S 
of p represents the single-qubit case. b(t) can be inter- 
preted as the amplitude of the upper level |e) of the qubit 
initially prepared with pf e (0) = 1 and 6(0) = 1, and is 
given by the inverse Laplace transform 



b(t) = J2f- 



■F(a) 
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where the parameter s is a complex number and F(s) — 
-^[/(t)] — Jo /(t) ex P( — st)dt with the correlation func- 



tion f{t) = Ei\9i\ e Ult = JcLoJ(Loy St and S {1) = 
ujQ—ujny The explicit form of b(t) depends on the specific 
spectral density of the reservoir 

As the density matrix should be of hermiticity, nor- 
malization, and semi-positivity, any pair of initial states 
could be written as 



pf(0) = 

(4(P) = 
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with \fi\ 2 < a(l - a), H 2 < m(1 - p), {0,v)€ C and 
< a (fi) < 1, (a, R. So we have 



> = , a\b(t)\ 2 ()b{t) 
r ' 1 /?*&*(£) l-a|6(£)| 2 
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P21 j ~ 1 ^*fo*(t) l-p|&(£)| 2 
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Using the definition of the trace distance, we obtain 



V s (pf (<),/£(«)) - |&(t)| V^C*)! 2 (« - /^) 2 + I/? - ^l 2 , 

In what follows, we consider the case that the qubit in- 
teracts resonantly with a reservoir with Lorentzian spec- 
tral distribution 
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with 70 the Markovian decay rate and T the spec- 
tral width of the coupling which has been widely 
employed in quantum optics [l^]. We may distin- 
guish the Markovian and the non-Markovian regimes 
using 70 and T: 70 < T/2 means the Markovian 
regime and 70 > T/2 corresponds to the non-Markovian 
regime. Therefore, using Eq. (4), we have b(t) — 
cxp(-R/2) [cosh(^/2) + (r/K)sinh(Kt/2)], 7o < T/2 



and 6(f) = exp(-Ti/2) [cos(k£/2) 
70 > T/2, with k 



+ (T/k) wx(Kt/2)], 
We may check that 
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in the non-Markovian regime, all local minima of \b(t)\ 
approach zeros at r™ ln = 2[nn — arctan(«;/r)]/K with 
n= 1,2,3, i.e., \b{r™ in )\ = 0. According to Eq. (7), 
all nontrivial trace distances own the same monotonic- 
ity, so 2? 5 (pf(C in ),pf(C in )) = 0. Therefore, Eq. (1) 
is reduced to M s = max V s (pf (r™ ax ) , pf (C™)) 

PI, 2(0) 

and the maximum taken over all pair of initial states is 
equivalent to finding a trace distance whose local max- 
ima are larger than those of others. In what follows, we 
adopt Vff (pf (t), pf (t)) for such a trace distance and 
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for the non-Markovianity of a single qubit case, where the 
summation is over all local maxima of T>fj- (pf (t), pf (*)) ■ 
We first prove the theorem below. 

Theorem: There exists a maximum trace distance \b(t)\ 
at any instant time in Eq. (7) when a = p = 1/2, |/3| = 
H = 1/2 and = 1. 

Proof: Suppose V s (pf (t), pf (*)) = \b{t)\d{t) with 



d(t) = ^/ |6(t)| (« - + |/3 - ^1 where d{t) can be 
taken as the distance between the points P\ (\b(t)\a,/3) 
and P2 (\b(t)\ p, v) , and Pi(x,y) (i = 1,2) denotes the 
points with x <E M and y E C. It is convenient to check 
that [\b(t)\ a - \b(t)\ /2] 2 /[|6(i)| /2] 2 + \(3\ 2 / (1/2) 2 - 1 = 



a(l — a) 



< 



0, which implies that the point 

Pi (\b(t)\ a, (3) or P 2 (\b{t)\ p,v) is in (or on the circum- 
ference of) the ellipse 



[x-|&(i)|/2] 2 /[|6WI/2] 2 + |y| 2 /(l/2) 5 
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with x € K and y G C. We know that the maximum 
distance between the two points in an ellipse is in between 
the two ends of the major axis. Since \b{t)\ < 1 (t > 0), 
d(t) reaches the maximum 1 only in the case of \b(t) \ a = 
\b(t)\ p = \b(t)\ /2, i.e., a = p = 1/2 and \(3\ = \v\ = 1/2 
and \P-v\ = 1. Consequently, V s (pf (*), pf (i)) will 
also reach the maximum |i>(£)|. I 

It is easy to check that any pair of initial states sat- 
isfying the conditions in above theorem definitely owns 
the same trace distance |6(t)|. Since the maximum of the 
trace at any instant time is for any pair of the ini- 

tial states not meeting the conditions in above theorem, 
the local maxima of \b(i)\ should be never larger than 
those of the initial pairs meeting the conditions. In what 
follows, we will employ Vff = \b(t)\ as the trace distance 
for measuring non-Markovianity. As a result, the calcu- 
lation of non-Markovianity can be simplified to an easily 
computable expression 
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with ax the time point when \b(t)\ reaches the nth local 
maximum. 



3 



(a) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



r 'T„ 

(b) (c) 




20 40 60 80 100 5 10 15 20 



V V 

FIG. 1: (Color online) (a) Blue bars (black in the printed 
version): The non-Markovianity of a single qubit coupled to 
a zero-temperature bosonic reservoir with a Lorentzian spec- 
trum. Light blue bars (gray in the printed version): The 
total growth of trace distance with the initial pair of states 
P| e )(0) = |e) (e| and p\ g )(0) = \g) (g\ . The maximum trace 
distance T>%- (blue solid lines) characterizes the decay pro- 
cess (red dashed lines) of a single qubit in (b) non-Markovian 
regime (e.g., T = O.I70) and (c) Markovian regime (e.g., 
r = IO70) respectively. 

Straightforwardly, we can find the relationship between 
the population of a single qubit initially in the excited 
state |e) and the maximum trace distance V\ e \ = (T>fj-) 2 1 
which yields 

= E \Av r ™ ax )' ( 12 ) 

n 

from which the non-Markovianity of the qubit coupled to 
the reservoir could be measured from the population of 
the upper level of the qubit. The requirements for this 
implementation are (1) the bosonic reservoir is initially 
in vacuum state and (2) the initial state of the qubit is 
prepared in the upper level. 

Figure 1(a) shows a comparison between the non- 
Markovianity (the dark blue bars) and P|g> (the light 
blue bars) in the case of non-Markovian regime ranging 
from T = O.I70 to T = 70, where p is measured 

by the trace distance V s (p\ e \(t), p\ g \(t)) = \b(t)\ 2 with 
the initial pair of states p\ e ){0) — |e) (e| (a = 1,(3 = 0) 
and p\ g )(0) = I.9) (g\ (/1 = 0, v = 0). Due to \b(t)\ < 1 and 

|6(i)| 2 < \b(t)\, the light blue bars are always shorter than 
the dark blue ones. In addition, it can be seen that the in- 
distinguishability of the states grows with the increase of 
reservoir bandwidth T. This might be interpreted as the 
non-Markovian character becoming less evident when the 
coupling between the qubit and the reservoir decreases. 

We have studied in Fig. 1(b) the population of the 
excited state V\ e ) and T>fj- in the non-Markovian regime 



(e.g., T = O.I70), which shows the population V\ e ) re- 
viving with the increase of T>^-. This could be explained 
as the non-Markovian effect of the bosonic reservoir: The 
process with T>fj- being larger corresponds to the case that 
information lost by the qubit flows back from the reser- 
voir, increasing the distinguishability. So the population 
revives for several times during this period. In contrast 
to the non-Markovian regime, no revival of population 
occurs in the weak coupling Markovian regime (e.g., V = 
IO70) as shown in Fig. 1(c), because T>ff monotonously 
approaches zero with b{i) ~ exp[— (r — K)t/2], This im- 
plies no information flowing back to the qubit. 

The above analysis can be conveniently extended to 
other spectral densities, other than Lorentzian spec- 
tral distribution, for non-Markovian characterization of 
a qubit coupled to a bosonic reservoir [l6| , as long as the 
condition |&(r™ in ) | = (n = 1, 2, 3, • • • ) is fulfilled [l7|. 

Let us take a brief look at the non-Markovianity of 
two identical non-interacting qubits A and B locally in- 
teracting with two independent zero-temperature bosonic 
reservoirs respectively. We noticed that taking the maxi- 
mum over any pair of initial states according to Eq. (1) is 
nearly intractable in two-qubit case, although numerical 
simulation might been employed for this job [Hj]. How- 
ever, since any growth of the trace distance is a clear 
illustration of non-Markovian character, in the following, 
instead of finding the maximum, we consider the trace 
distance 

V T (p l++) (t),p h _ } (t)) = \b(t)W2-2\b(t)\* + \b(t)\* 

as the measurement of the non-Markovian character of 
the two qubits, where the initial pair of states P\++) (0) = 

\+)a (+1 ® \+)b (+1 and P\— >(°) = \~)a (-1 ® \-)b (-1 
with |±) i = (|. 9 ). ± |e) i )/V5 and i = A,B. Similar to 
the single-qubit case, the trace distance approaches zeros, 
i.e., V T {p ]++) {TT%P\—){rt n )) = when \b(r^)\ = 
0. Therefore, a lower bound of non-Markovianity can be 
calculated by 

< ++ >,pi-> ^ E l fo (C aX )l v/2 - 2|6(r„- x )| 2 + \b{r^)\\ 

n 

(14) 

or by the equivalent form 

K ]++) ,pi- } = E V 27? i«>( T " ax ) - 2 ^>( r " iax ) 2 +^ie)(r I r x ) 3 - 

n 

Similar to the single-qubit case, the non-Markovian char- 
acter lowers with the increase of spectral width of the 
coupling, as shown in Fig. 2(a). 

Conventionally, the dynamics of two entangled qubits 
can be measured by Wootters concurrence [l8|]. Using 
the method in Ref. [l{|, we obtain the concurrence of 
the two qubits as 

C m = max{0,|fe(t)| 2 }, 

C m - max{0,|6(i)| 4 }, (16) 
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FIG. 2: (Color online) (a) A lower bound of non-Markovianity 
corresponding to the trace distance T> T with initial pair 
of states P]++) (0) = \+} A {+\ <S> |+) fl (+| and p, — > (0) = 
\-)a(-\®\-)b(-\ (l±)i = (\9) i ±\e) i )/V2,i = A,B). The 
trace distance T> T (blue solid lines) characterizes the dynam- 
ics of two non-interacting qubits coupling to two indepen- 
dent reservoirs in (b) non-Markovian regime (e.g., F = O.I70) 
and (c) Markovian regime (e.g., V = IO70). While the 
red dot-dashed lines and green dotted lines represent the 
concurrence of two qubits initially prepared in Bell states 
I*) = (|ffe> + \eg))/V2 and |$) = (\gg) + \ee))/y/2, respec- 
tively. 



when the initial states are prepared in Bell states I*) = 
jj|e) + \eg))/V2 and |$) = (\gg) + |ee))/\/2, respectively 

The trace distance V T (Eq. (13)) measuring a lower 
bound for non-Markovianity can well characterize the 
non-Markovian behavior of the dynamics of two qubits. 
In Fig. 2(b), the dynamics of the qubits initially prepared 



in Bell states |\T/) (red dot-dashed line) and |<I>) (green 
dotted line), respectively, in non-Markovian regime (e.g., 
r = O.I70) is investigated. The entanglement measured 
by concurrence periodically vanishes in accordance with 
the trace distance T> T . In contrast, the trace distance 
T> T , in the Markovian regime (e.g., T = IO70), asymp- 
totically approaches zero, as shown in Fig. 2(c) which 
implies that no entanglement revival exists in Markovian 
regime. 

In practice, our scheme would be very preferable for ex- 
perimental implementation. According to Eqs. (12) and 
(15), we can study the non-Markovian effect on coher- 
ence and entanglement of the qubits coupled to bosonic 
reservoirs by directly measuring the population of a qubit 
without resorting to tomographic reconstruction of the 
density matrix. Besides, our proposal requires no spe- 
cific information about the interaction between the qubit 
and the reservoir. These favorable features make our 
scheme feasible in experiments under real environments, 
e.g., using two-level atoms confined in optical microcav- 
ities [20[ or under simulated reservoirs [2l|, e.g., using a 
spin- reservoir model with spectral densities ranging from 
sub ohmic to super ohmic cases simulated by trapped 
ions HI]. 

To summarize, we have presented a simple method 
for measuring the non-Markovian character of the qubits 
coupled to bosonic reservoirs. We believe that this eas- 
ily operated measure for non-Markovianity would be very 
useful for further understanding non-Markovian behavior 
and also for experimental exploration, particularly in a 
realistic experimental situations without knowing much 
about the interaction between a qubit and the environ- 
ment. 
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